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showed higher amplitudes than expected on the basis of their t* data. Still, 
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This report presents the results of the second phase of the Nevada Test 
While stations in the Western United States (WUS) are 
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E Site (NTS) experiment. 
a no longer confined to the Nevada Test Site (NTS) in this phase of the work, 


E the name "NTIS experiment” will remain in use for the project, even though a 


more proper name might be “reciprocal attenuation measurement project." In 


this phase the recording stations were moved to locations shown in Figure 1. 


3 At NTS a second station, OB3NV, was set up next to OB2NV. OB2NV remained in 


its previous position on the Climax Stock, and the Pahute Mesa stations were 
moved to the Yucca Flats, a broad alluvial valley south of the Climax Stock. 


Four sensors (YFNV, YF2NV, YF3NV and YF4NV) were set up in the Yucca valley in 
In 


Ee ee 


positions shown in Figure 2; all these stations were recorded digitally. 
additien to the NTS stations, two digitally recording stations, RKON and HNME, 
To extend the number of sites sampled, Special Data Collec- 


nea 


remained in place. 
tion System (SDCS) stations were set up at the sites of the FAULTLESS explosion 


SS Sn 


(FANV) in Nevada and the GASBUGGY explosion (GBNM) in New Mexico, and also 


above the Tatum salt dome in Mississippi (TQMS), the site of the SALMON nuclear 
The approach 
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5 explosion. FANV, GBNM and TQMS recorded on analog magnetic tape. 
and methods used in data analysis are similar to those used in the first phase 
of the NTS experiment described in Seismic Data Analysis Center (SDAC) TR-77-7 


‘ graphical coordinates of all the stations 
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and TR-77-9 (Der et al., 1977a,b). 
discussed in this report are given in Table VII. 
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: In addition to discussing data obtained from the second configuration of 


‘igs HYP 


seismic stations, we will interpret it, along with data from the first con- 


figuration; the same will be done with the total data set when the project is 


completed. The experiment provides individual dm, and At* values for a fairly 


Der, Z. A., M. S. Dawkins, T. W. McElfresh, J. H. Goncz, E. G. LaPella, and 
M. D. Gillispie, 1977a, Teleseismic P wave amplitudes and spectra at NTS | 
and selected Basin and Range sites as compared to those observed in 
Eastern North America, NTS experiment - Phase 1, Final Report; SGAC-TR- 

71-7, Teledy.e Geotech, Alexanaria, Virginia. 


Z. Aw, M. S. Dawkins, T. W. McElfresh, J. H. Goncz, E. G. LaPella, and 
M.D. Gi*Lispie, 1977b, Teleseismic P wave amplitudes and spectra at NTS 
and SHOAL site as compared to those observed in eastern North America, 
Preliminary Report; SDAC-TR-77-9, Teledyne Geotech, Alexandria, Virginia. 
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large set of stations. These, and the estimates of crustal amplification 


Ay» facilitate study of inter-relationships between Am, At*® and Ay: 
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DATA ANALYSIS 


General 
With few refinements, analytical tools used in Phase II closely fvllow 


those used in Phase I (Der et al., 1977). For example, magnitude residuals 
were determined by a least squares adjustment procedure over the SDCS station 
network instead of taking direct differentials. The results of the least 
squares scheme are more reliable because all data were taken into account 
rather than daca limited to the common readings at the station pair under 
study, but the results are modified only slightly by the least squares pro- 
cedure. In addition to the standard Haskell matrix method for computing the 
crustal response, finite difference (FD) calculations were also used to esti- 
mate the effect of geological structures at Yucca Flats. (The added refinements 
do not affect the basic conclusions of the report.) Part of the analyses, the 
ca}. ulations of t* for short period S waves, are included in a separate report 
(Smart et al., 1978). Events used for analysis, along with the amplitude 


and period readings and the computed magnitudes, are listed in Appendix A. 


Epicentral data were taken from various lists such as the Network Event 
Processor (NEP) bulletin and the Hagfors array bulletin. This was necessary 
because searching the data (especially digitally recorded data) was impractical 
because of limited accessibility to the PDP-15 computer and the need to have 
some epicenters determined and available for data reduction. The epicentral 
distance range was limited to 25° < A < 85°, thus eliminating the effect of 
complications from upper-mantle travel time triplications at smaller distances 
and the effects of the core mantle boundary beyond A = 85°. “Film playbacks 
of the analog stations were used extensively for amplitude readings because 
of the limited capability of analog to digital conversion on the PDP-15,. 


Selected signals were subsequently digitized for spectral analysis. 


Smart, E. Z. A. Der and A. H. Chaplin, 1978, Short period S wave attenuation 
under North America, SDAC-TR-78-6, Teledyne Geotech, Alexandria, Virginia, 
in preparation. 
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Body Wave Magnitude Measurements 


j 
| 
{ 
In past SDCS experiment reports, biases between stations were computed | 
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directly using the magnitude differential average from selected key station 
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pairs. While this approach was less complex, it did not take into account 


measurements that involved only one station of the pair in question. However, 


even though operational difficulties (non-overlapping station life~times) 


imposed epicentral distance criteria, and highly variable noise conditions 


hindered P amplitude measurements for most events at all stations, all measure- 
a ments should optimally be utilized. j 
: ; 
Be 2 
g m,, differentials between the stations was found to be a strong function of 


distance between stations Physically, this stemmed from loss of signal 


: coherence with relative distance resulting from multipathing, source radiation 


patterns, and changes in crustal structure. Figure 3 shows a plot of the 


r standard deviation of a single Mm differential measurement versus the logarithm 


q of epicentral distance (in degrees); the relationship appears to be linear. 


5 Table I is a tabulation of data points in this figure. The total variation 


F 
E over the range of relative distances between stations is five-fold, which is 


During the first phase of the project, the standard deviation of relative | 
quite large. Thus, the accuracy of An, will largely depend upon distances 
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between stations, 
: To incorporate all measurements into determining magnitude differentials 


4 across the SDCS network, as well as to take into account the effect of signal 


coherence as manifested in the Ohi - A° relationship, a least squares scheme 
has been devised. For a given event we write the magnitude differentials 


between pairs of stations i and j in the form 


ie hh are 
A 


i j ._k, ij o 
Am, - Am = Am, + €,,0 ij) (1) 
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where is the station term as fe) station 1, is the observe 
am, i \ i (bias) of i i, “om fe wheeeved 
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magnitude differential between station i and j for event k and €,,(A°4j). an 
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error term dependent upon the distance A°ij between the stations. 


pected value of this error term 
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TABLE I 


Standard deviation of the magnitude differentials versus 
interstation distance. 


Am, 
° Amp least single 
Spatton att a "amp N squares diff. diff. 
OB2 - OB3 0.00789 -095 113 ~0.016 -0.016 


YF - YF2 0.00768 - 069 47 ~0.064 ~0.077 
YF - YF3 0.0158 - 104 33 -0.128 ~0.095 
YF - YF4 0.0318 -091 36 ~-0.107 -0.079 
YF4 ~ OB2 0.161 ~283 35 0.413 0.413 
GB — OB2 7.09 ~357 89 0.058 0.062 
FA ~ OB2 1.41 -319 108 -0.020 -0.014 
FA - GB 7.39 432 160 -0.078 -0.063 


HN - RK 17.61 .385 71 -0.040 -0.099 
RK — OB2 20.97 456 167 0.156 0.146 
HN - OB2 36.50 505 58 0.116 0.141 
YF - FA 1.58 353 42 0.326 0.321 
YF - GB 7.04 400 30 0.248 0.208 
RK - FA 20.04 401 52 0.176 0.150 
RK - GB 17.16 .530 46 0.098 0.126 
RK ~ YF 21.05 446 32 -0.150 ~0.223 


YF4 - FA 1.57 - 268 31 0.433 0.440 
YF - OB2 0.17 - 303 50 0.306 0.304 


0 = .310 + .107 10g, 4° 


om, 


Goodness = .9 


Np = 18 
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can be read from the regression line in Figure 3 as a function of distance. 


Because taking differences in all possible combinations is redundant, a 
hierarchy of stations was assigned, and we used the leftmost available station 
bias term as the reference (positive) term in equation (1). The hierarchy 
chosen was 


0B2, OB3, YF, YF4, YF2, YF3, NT, NT2, FA, GB, RK, HN. 


It was selected to optimize the distance distribution so that most of the 


distances used, and consequently values of 0, » are as small as possible. 
Thus OB2, at the center of the NTS cluster and Yocated between GB and FA, is 


the prime candidate for the most commonly used reference station, while RK 
and HN are outliers and rank low. We also chose Am, = 0 for OB2, which merely 
defined the reference for other Am, (we could have chosen any other station). 
Magnitudes computed without division by the dominant period T; (called m) 


were adjusted in a similar maaner. 


The results are summarized in Figure 4, and in Table II, which consider 
all data from the first two phases of the NTS experiment and also compare the 
SHOAL site (SZNV) with the Sleepy Eye, Minnesota (SEMN) site. The SZNV-SEMN 
comparison was placed within the framework of the NTIS experiment by assuming 
that RKON and SEMN possess identical atteruation and amplitude properties. 

This assumption was reasonable because both sites are located on the Canadian 
Shield granite. All quantities in Figure 4 are relative to the station OB2NV. 
The top of the figure shows the raw magnitude residuals. At stations located 

on the shield (HNME, RKON and SEMN), m, values seem to be high relative to OB2NV. 
At NTS, stations on Yucca Flat and the Pahute Mesa are also high, while RANV, 
GBNM, and OB3NV are at about the same mm, level as OB2NV. SZNV (assuming the 
equivalency of the SEMN-RKON pair) is low in mr relative to OB2NV. The plot of 
quantities ms (magnitude without division by T.) follows the pattern of Am, 
indicating that amplitudes primarily determine Am, variation despite the fact 


that the RKON-OB2NV differential in ms becomes statistically insignificant if 


b 


Directly determined m mi and T. differentials for selected stations ate 


the division by T, is removed. The values of Am’ are not available for SZNV. 


also shown, as histograms, in Appendix B. Because these are not adjusted in 
the least squares sense, they are somewhat different from numbers given in 
Table IL. 
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Figure 4. 
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TABLE IL 


Magnitude differentials Am, Am‘> AT, and om Ort and t* relative 


to OBZNV for the stations discussed in this report. 


FA -.020 + .061 -.055 + .074 -.036 + .054 ~.020 + .068 -.288 


Am, Am’ At* AT 1 Sai 

WME = «126 + 6119 = 100 + 1440-078 + .051 NA. 116 

SEMN : 

RKON «156 + .071 «055 + .086 «= -.187 + 015 = -.192 + .038 =. 156 

YF 1306+ .056 «294 + .068 = -.040 + .050 ~.051 + .071 046 | 

YF2 «370 + 065.367 + 079 «O14 + 055 =. 011 + 1067-050 | 

YF3 «434+ ..076 «= .43€ + 092.016 + .063 032 + 074-164 | 

YF 413 + .080 403 + .097 «006 + 062. ~.013 + .073 «113 | 
| 
| 


GB 058 + .074 .033 + .089 .013 + .050  -.068 + .052. -.252 
NT 246 + .072 186 + .087 -.070 + .060  -.100 + .060 ~—. 03 
NT2 .237 + .057 175 + .069 = -.060 + .050-- -.100 +. 05 027 
OB3 .016 + .018 .016 + .023 -.036 + .019 .007 + .027. 016 

S7 ~.064 + .09 N.A. -.017 + .096 NLA. -.064 

; 
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Possible Bias Effects in m Due To Variable Thresholds and Noise Levels 
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The criteria set for amplitude measurement required that S/N be greater 


than 3. Thus, event pairs were eliminated when one of the stations had S/N < 3, 


Because these measurements could be those with the low signal at station A 
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relative to station B, the average amplitude at A could be overestimated 
relative to B at low magnitudes (Herrin and Tucker 1972, von Seggern and 


Blandford 1976). Only if the great majority of all readings is considerably 
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4 


above the 3:1 S/N level can such bias in the average be ruled out. We have 


chosen to plot the differentials of m, for selected station pairs against the 


SSN f 


3 average mi, for the same two stations. Pronounced trends in plots such as these 
a indicate biases in procedures used for determining Am, . A slight change of 


| 4m, | with increasing magnitude might also suggest a shift of corner frequency of 
This shift would be more visible at a : 


seismic sources to lower frequencies. 
Absence of a clear trend indicates that bias 


wea eae 


high Q than at a low Q station. 
in Am, from variable noise levels is not significant. 

Figures 5 through 11 show such plots for a selected set of key station 
None of these plots, including the most critical pair RKON-OB2NV, 


auc 


pairs. 
suggests a clear trend and, therefore, bias effects from our procedure are 


probably negligible. Note that the noise level is approximately 2 times higher 
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at RKON than at OB2NV. However, the raw amplitudes on the film (magnitude at 


AS 


1 Hz, not amplitude corrected at T, or A/T) also average 2 times higher at RKON. 
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Thus, the average expected S/N is the same and no bias would be expected. 
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Herrin, E. and W. Tucker, 1972, On estimation of body wave magnitudes, Report 
to the AFOSR: Dallas Geophysical Observatory, Southern Methodist Univer-— 


sity, Dallas, Texas. 
von Seggern, D. H. and R. Blandford, 1976, Seismic threshold determination; 
Bull, Seism, Soc. Am., 66, 753-788. 
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Figure 5. Am, versus averaged magnitude for the station pair RKON-OB2NV. 
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Figure 6. Am, versus averaged magnitude for the station pair HNME-OB2NV. 
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Am, versus averaged magnitude for the station pair FANV-OB2NV. 
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Figure 8. Amp versus averaged magnitude for the station pair GBNM-OB2NV. 
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Figure 9. Amp, versus averaged magnitude for the station pair NT2NV-OB2NV. 
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Figure 10. Amp versus averaged magnitude for the station pair NTNV-OB2NV. 
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CALCULATIONS OF THE CRUSTAL RESPONSE 


In this report, as in the previous NTS experiment report, we attempted 
to estimate the crustal contribution to the observed magnitude bias at each 


observing station. The Yucca Flats sites rest on thick unconsolidated sedi- 
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ments and tuff that cause considerable signal amplification (Houser, 1968; 
Fernald et al., 1968; Healy, 1968; Ramspott and Howard, 1975; Hays and Murphy, 
1971). The FANV site is also located over alluvium and tuff, but according 

to test site information, the alluvium is more consolidated at this site 
(McKeown and Dickey, 1969); (Lt. Col. George Bulin of ARPA also provided us 


with data relevant to the FANV site). Alluvium and a thick sedimentary car- 


bonaceous~shale sequence also underlay the GASBUGGY site (Thornbrough, 1971). 


Because the part of the structure that primarily determines crustal amplifi- 


cation was found near the surface, the structures were modeled only down to 


the basement (Der, McElfresh and Mrazex, 1977), which was assumed to have the 


same elastic properties as the granite stock at OB2NV that was modeled by a 


simple homogeneous halfspace. This similarity makes the pulse sizes of the 
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Houser, F. N., 1968, Application of geology to underground nuclear testing, 
Nevada Test Site, Geol. Soc. Am. Memoir, #110, E. B. Eckel, Editor, 
Boulder, Colorado. 


Fernald, A. T., G. S. Corchary, W. P. Williams and R. B. Cotton, 1968, Sur- 
ficial deposits of Yucca Flat area, Nevada Test Site, Geol. Soc. Am. 
Memoir, #110, E.B. Eckel, Editor, Boulder, Colorado. 
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Healey, D. L., 1968, Application of gravity data to geological problems at 
Nevada Test Site; Geol. Soc. Am. Memoir, #110, E. B. Eckel, Editor, 
Boulder, Colorado. 


Ramspott, L. D. and N. W. Howard, 1975, Average properties of nuclear test 
areas and media at th USERDA Nevada Test Site, Lawrence Livermore 
Laboratory, -UCRL-51948. 
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Hays, W. W. and J. R. Murphy, 1971, The effect of Yucca fault on seismic wave 
propagation; Bull. Seism. Soc. Am., 61, 697-706. 


McKeown, F. A. and D. D. Dickey, 1969, Fault displacements and motion related 
to nuclear explosions; Bull. Seism. Soc. Ams, 59, 2253-2271. 


Thornbrough, A. D., 1971, Current practices and anticipated improvements in 
PNE, In "Peaceful Explosions II", by the International Atomic Energy 
Agency, Vienna. 
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computed synthetic records comparable without correction for halfspace pro- 


wy) 


perties below the layered stock. Pulse sizes of a 50 kt explosion, as modeled 
by von Seggern and Blandford (1972), were compared after passing them through 


our layered halfspace models. We also attenuated each pulse with a multi- 


Hits a 


plicative spectral factor exp (-wft*) where t* was chosen to be .45, a typical 


Tea 


value for the WUS. By removing most of the high frequencies, the attenuation 


Nene Ree 


rT 


; factor makes the pulse more rounded. This pulse has a spectrum which, in 

4 spectral content, well represents the average teleseismic P-wave arrivals. 

Ee It is peaked at 1 cps. and falls off at a rate of somewhat more than w “ at 

a high frequencies. To reduce variations caused by changes in the angle of in- 
cidence, we computed synthetics for three angles (20°, 25°, and 30° measured 

= from the vertical) and then averaged the relative amplification factors between 


stations obtained for these three angles. 


é Table III lists the model parameters used in the calculations, and Figures 
a 12 through 18 detail the results of the calculations. All the figures are 

laid out identically; the model name is at the top, and the angles of incidence 
are on the left followed by the crustal amplitude responses plotted as functions 
of frequency. In the center are the impulse responses of the crustal model, 

and at the right are the synthetic seismograms plotted on the given time scale, 


All plots are normalized to unity in maximum, and the numerical value of the 


normalization factor is indicated on the right of each plot. The values of 


E crustal amplification are summarized in Table IV. 
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von Seggern, D. H. and R. R. Blandford, 1972, Source time functions and spectra 
for underground nuclear explosions; Geophys. J. R. Astr. Soc., 31, 83-97. 
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TABLE IIT 


Crustal models used in the computation of correction 
factors by the Haskell matrix method. 


GBNM 
d o g § 
700 2.00 1.06 2.00 
2150 3.25 1.63 2.00 
.200 3.45 1.73 2.00 
.120 3.90 2.16 2.28 
-600 4.40 2.48 2.49 
625 4.80 2.73 2.55 
© 5.70 3.29 2.70 
FANV (McKeown and Dickey, 1969) 
0.44 — 2.5 1.06 2.3 
1.00 3.0 1.60 2.3 
1.00 3.5 1.73 2.3 
1.525 4.0 2.20 2.7 
Yucca Flats (Hays & Murphy, 1971) 
d d d d 
YENV YF2NV YF3NV YF4NV a B 6 
18 24 29 29 1.3 .659 1.75 
295 58 61 70 2.0 1.07 1.196 
oo cod cod fad 5.7 3.36 2.7 
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In addition to calculating cru..tal response using Haskell's matrix method 
for horizontally layered media, a finite difference method for inhomogeneous 
media was also used to estimate crustal amplification at the four stations on 
Yucca Flats (the method is described in the paper by Kelly et al., 1976). We 
atten; “ed to model Yucca Flats with a structure derived from tiays and Murphy 


(1971), utilizing Ramspott's and Howard's (1975) and Fernald et al.'s velocity 


and structural data. 


Figures 20a-d together show the vertical displacement wavefield at various 
times caused by a single cycle P wave incident on the Yucca Flats model. (The 
method for making this wavefield js outlined in Figure 19.) Assume that an in- 
cident wave with displacement amplitude A, plus signs, denotes upward displacement 
with amplitude [A| > .33; and that minus signs denote downward displacement with 
similar amplitude. The incident wave has a frequency of 1 Hz. The figure also 
indicates geological units, somewhat oversimplified, and the location of Yucca 
fauLt and the four stations at Yucca Flats. To increase computational sta- 
bility, transition zones were inserted between the various geological forma- 
tions. The "snapshots" of vertical displacement wave fields in Figures 20a, 

b, and c show, in successive order, the slowing down of the wavefront by the 
basin (20a) and the development ox the in-phase reflection from the surface 
(20b and c). Figure 20d shows the synthetic wavetrains at five stations, four 
at Yucca Flats and the fifth at an imagin.ry reference station on Paleozoic at 
the side of the valley. The figure shows that the time delay at all Yucca 
stations is the greatest at YF4NV where the sediments and volcanics are the 
thickest. The amplification of the waves at Yucca Flats is also evident re- 
lative to the reference station. Since we needed to compare Yucca Flats to 
OB2NV, a station on granite, we used a granite halfspace for a reference model 
in computing the amplification at Yucca Flats, instead of the reference station 
in Figures 20a-d. The base 10 logarithms of the amplitude ratios were used 


for the crustal magnitude corrections; they are summarized in Table IV. 


a 
Kelly, E. R., R. W. Ward, S,. Treitel and R. M. Alford, 1976, Synthetic seis— 
mograms, 2 finite difference approach; Geophysics, 41, 2-27. 
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TABLE IV 


Summary of crustal amplification results. 


Stations Magnitude differential relative to OB2NV 
Finite Difference layered Am, 
1Hz (Haskell) observed 
YF -26 34 .306 
YF2 32 34 ~370 
YF3 .27 232 434 
YF4 . 30 232 2413 
-39- : 
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Explanatton of the method used in Figure 20. 


Figure 19. 
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Figure 20d. Finite difference calculations for Yucca Flats. 
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siderable amplification at Yucca Flat, and the two methods yield similar esti- 


Both the Haskell matrix and the finite difference methods suggest con- 
| 
H 


mates. Note the reverberations of the waveforms at Yucca Flats in the synthetics 
computed by the Haskell method. The lengthening of the P wavetrains due to : 


reverberations can also be seen by inspecting the recorded actual waveforms. 
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SPECTRAL ANALYSIS 


Relative t* values for the SDCS stations have been computed. Because com 
puting spectral ratios for all possible pairs of stations is redundant, we 
computed ratios only for those station pairs directly connected by lines in 
Figure 21. Differentials in t* (At*) obtained from spectral ratios and their 
standard deviation for other station pairs can be easily derived from these 
figures. This approach takes advantage of the fact that, for closely located 
station pairs, aver has smaller variance than for distant pairs of stations. 

For example, determining At* for the RKON-OB2NV and OB2NV-YFNV pairs is more re- 
liable than determining At* directly from the fewer events common to RKON and 
YFNV alone. This statement is true because many events for the RKON-OB2NV pair 
were available during the project's two phases to reduce the variance of their 
mean At*, while few events were sufficient to define At* for the OB2NV-YFNV pair. 
The histograms of the measured At* are given in Appendix C. In addition to 
those involving the new stations of Phase II, updated versions of histograms for 
the OB2NV-RKON and HNME-RKON pairs are presented. In all cases, an accuracy 

(26 4 4) of ~ .05 sec was the goal. 


The easiest way to discuss At* is to compare them to a common standard 
station. Some of these At* mean values and their standard deviations were derived 
indirectly. The summary figure (Figure 4) uses OB2NV as a standard station. 

The figure shows that all WUS stations have 2ssentially the same t* as OB2NV with 
the exception of stations NINV and NT2NV, which have significantly lower t* than 
OB2NV. The OB3NV-OB2NV differentiai is also significant statistically but the 
difference in t* so small that it is of no importance in this study. The averages 
of the differential in dominant period (AT) show variations similar to At*, in- 


dicating that variations in frequency content are visible on the time domain 


traces. The RKON-OB2NV differential in t* is about .2 sec and highly significant 


statistically. The HNME-OB2NV differential is less (.08 sec), but it is also 


significant at the 95% confidence level. This lower value may indicate some 
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attenuation under the Northeastern U. S.; a possibility suggested in Solomon's 


and Toksoz's (1970) work and Der et az.'s work (1975). 


The Yucca Flat stations had some difficulty in reliably determining 
t*, This difficulty resulted from high, unstationary man-made noise (drilling) 
at these sites that made the standard procedure for estimating noise in the 
signal window (taking spectra of time windows prior to the arrival of the P 
wave) fail repeatedly. Many of the spectral ratios supposedly show points 
with high S/N (around 4 Hz) that disagree with the trend at lower frequencies. 
This can be attributed to nonstationary noise rather than to real signal energy 
at 4 Hz, which should be far beyond the corner frequency expected for most 


events. Crustal response calculations rule out major enhancements of high 


frequency energy. (Thus, the reader should assign less weight to the t* deter- 


mination at YF stations). 


The standard deviation of relative t* values depends upon mutual distance 


between stations similar to that shown by the magnitude residuals. Table V 


includes a tabulation of o versus A° as shown in Figure 22. Note that the 


At* 
numerical values of Ones are smaller and increase more slowly with distance 


| 


than on - Because absolute numerical values of Am, and At* are both of the 


order of a few tenths in magnitude units and seconds, respectively, fewer 


measurements are needed to establish At* between stations than to determine 


alk te 


magnitude residuals with the same numerical accuracy in the respective units 


used in this report. This result is in agreement with experience at seismic 


Mg 


Tenor 


arrays and shows that At* is numerically a more stable quantity than 4m, in 


terms of multipathing and other disturbances (Der et al., 1977b). 


Solomon, S. C. and M. N. Toksoz, 1970, Lateral variations of attenuation of 
P and S waves beneath the United States; Bull. Seism. Soc. Am., 60, 
819-838. 


Der, Z. A., R. P. Masse and J. P. Gurski, 1975, Regional attenuation of short— 
period P and S waves in the United States; Geophys. J. R. Astr. Soc., 
40, 84-106. 

Der, Z. A., M. S. Dawkins, T. W. McElfresh, J. H. Goncz, E. G. LaPella, and 
M. D. Gillispie, 1977b, Teleseismic P wave amplitudes and spectra at NTS 
and selected Basin and Range sites as compared to those observed in 
eastern North America, NTS experiment - Phase I, Final Report, SDAC-TR- 
77-7, Teledyne Geotech, Alexandria, Virginia. 
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TABLE V 


Standard deviation of At* versus interstation distance. 


Station Pair a “ate x 

GB - OB2 7.09 2222 25 

FA —- OB2 1.41 -139 27 

YF4 - YF .0318 -098 29 

YF4 — YF2 -924 - 066 25 

YF4 - YF3 26558 .078 26 

YF3 - YF - 0138 .091 26 

YF3 - YF2 - 60768 -093 22 

YF2 - YF - 00768 - 069 36 

YF —- OB2 17 -155 39 

P OB3 - OB2 - 00789 - 065 4Eé 
z OB2 - RK 20.97 -149 75 
: RK —- HN 17.61 -140 24 
NT2 - NT -100 2125 27 

NT2 -— OB2 -192 -156 38 

NT2 — RK 21.03 -212 43 

Gaps = 0-147 + .032 log, 96° 
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Joint Consideration of 4m, An’, At*® and the Crustal Corrections 


: corr F af : 
The bottom of Figure 4 shows the 4a, magnitude residuals relative to 


OBZNV, corrected for estimated crustal response; the stations on granite or 
hard metamorphic rock are left uncorrected (HN, RK, SE, SZ, 0B3). Stations 
were correct.d on the thick tuff sequence of Pahute Mesa (NT, NT2), the allu- 
vium and tuff of Yucca valley (YF, YF2, YF3, YF4) and sediments of other kinds 
(FA, and GB). Ideally, a negative cvrrelation should exist betwe2n a 
and t*, if t* and the crustal corrections are the sole factors determining 
m- however, tnis correlation does not hold because, while estimated crustal 


corrections reduce most NTS stations to roughly the same magnitude level (with 


the exception of YF3NV and YF4NV), FANV and GBNM when corrected fall below the 


NTS levei. Near surface information is quite reliable at the FANV site, where 
several boreholes have been logged for velocity, but the accuracy of GBNM 
corrections is quescionable. Still, available information suggests that re- 


ducing the FANV and GBNM voints below the a level of NTS appears unavoidable. 


Apparently, then, the consistency between the set of values of Am,» Amt > 
At* and ‘he crustal corrections for stations in Phase I was destroyed after 


adding sew data. The conseyeences of this finding are discussed below. 
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Interpretation of the Results 


The purpose of the NIS experiment is to measure the amount of anelastic 
attenuation in the short—period band under a selected set of nuclear test 
sites within the continental United States by measuring spectra and amplitudes 
of incoming P waves from teleseisms. Regional patterns of magnitude diff- 
erential measurements for teleseisms have been previously vutlined (Guyton, 


1966; Evernden and Clark, 1970; Booth, Marshall and Young, 1974; Der, Massé 


Evernden, J. and D. M. Clark, 1970, Study of teleseismic P. II amplitude data; 
Phys. Earth Planet, Interiors; 4, 24-3. 


Booth, D. C., P. D. Marshall, and J. B. Young, 1974, Long-and short~period 
amplitudes from earthquakes in the range 6° ~ 114°; Geophys, J. R. Astr. 
Soc. 
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and wturski, 1975; North, 1977), and little doubt can exist about the presence 
of regional differences in m, measurements. While not a central part of the 
experiment, we observed regional differences in the frequency content of P 
waves that appear to correlate with the magnitude residuals (Der and McElfresh, 
1977). On an average regional level, the magnitude spectral data and crustal 
amplification estimates also seemed to correlate (Der, McElfresh and Mrazek, 


1977) and a similar regional pattern seemed to exist for short period S waves 
(Smart, Der and Chaplin, 1978). 


wisn 


While on a regional basis the correlation between At*, Am» and crustal 


ry 


corrections appears statistically significant, determining the degree of an-~ 


elastic attenuation under a specific location is quite different. Therefore, 


ay 
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‘ 
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although critical measurements for estimation of anelastic attenuation at a 


3. 
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E given site include spectral slopes and body waves amplitudes, exactly which 

Q measurement more effectively diagnoses the reciprocal Q effect of the mantle 

3g under the given site should be decided by careful evaluation of the reliability 

g of the diagnostics used. Both spectra and amplitudes are subject to effects 

a unrelated to anclastic Q, and we want to base our calculations on the most 

q stable and diagnostic combinations of measured quantities available. Therefore, 

a the choice of est_mators for upper mantle Q should be based both on experience 
4 and theory. 


Guyton, J. W., 1964, Systematic deviations of magnitude from body waves at 
seismograph stations in the United States, Proc. VESIAC Conf. Seismic 
Event Magnitude Determination, University of Michigan, 4410-71-X. 


3 


North, R. G., 1977, Station magnitude biases~its determination, causes and 
E. effect, Lincoln Laboratory, M.I.T. Technical Note 1977-24, 62. 
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Der, Z. A. and T. W. McElfresh, 1977, The relationship between anelastic 
attenuation and regional amplitude ano- ilies of short period P waves 
in North America; Bull. Seism. Soc. Am., 67, 1303-1317. 


a Der, Z. A., T. W. McElfresh and C. P. Mrazek, 1977, The effect of crustal i 
3 structure on the station magnitude anomalies (magnitude bias), SDLAC- j 
TR-77-1, Teledyne Geotech, Alexandria, Virginia. 
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Smart, E., Der, Z. A. and A. H. Chaplin, 1978, Short period S wave attenuation 
under the United States, SDAC-TR-78-6, Teledyne Geotech, Alexandria, 
Virginia. 
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In the absence of other effects, slopes from the ratios of body wave 
spectra are useful in evaluating anelastic attenuation. The attenuation cor- 


rection for Mm, should be, in the first approximation, proportional to t* 


pe (1) 


Numerical simulations can provide more exact relationships. Inclusion of the 
effect due to the local crustal amplification factor Aw leads to a more real- 


istic tormula for m, at a given station 


™, z a+b logs oA, +ect*® +e (2) 


where a, b, and c are constants. a is a simple additive term, b should ideally 
be unity if A. is unbiased, and c could ve derived from numerical simulation. 
Finally, € is an error term that should be zero, if all other effects of mag~ 


eas 


nitudes are negligible, and if A. and t* are accurate. If e is large, cannot 
be used to measure attenuation. Of all quantities in equation (2), t* is most 


important in terms of determining magnitude bias for explosions at teleseismic 
recording stations. If known, then this bias determines the reciprocal upper 
mantle attenuation. The question is whether measuring m, at individual stations 


improves the definition of t*. 


Experience indicates, however, that amplitude measurements are almost use- 
less in determining t* at an individual station. This is true because of our 
crude knowledge of crustal structure at most locations. In addition, horizontally 
homogeneous models used for computing A. are not realistic, and therefore, A, 


may have a large error associated with it. 
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The teleseismic magnitude residuals (Booth, Marshall and Young, 1974) 


at selected stations of the LRSM network provide a set of empirical data 


Suitable for determining the coefficients in equation (2). Using available 


geological information and the LRSM site reports, horizontally layered crustal 


models were constructed for each site. The crustal amplification was 


estimated by passing the normalized pulse of a 50 kt explosion through the 


layered stack and then, with Haskell's (1962) algorithm, computing the sur- 


face amplitude. The Am,, logarithms of crustal amplification factors, and 


b? 
t*, (.2 for EUS, and .45 for WUS) based on our study of regional differences 


(Der and McElfresh, 1977), were fitted with the linear model of equation (2). 


(Details of the analysis were given in a previous report (Der, McElfresh 


and Mrazek, 1977), so only the major conclusions are repeated here.) 


Regression analysis of the data shows that both "b" and "c" are statistically 


significant with the values and 95% confidence limits given as 


b= 1.03 + .43 
ec = -1.35 + .32 


il 


The crustal effect for outgoing waves should be considered separately. 


Figure 23 shows a plot of the magnitude residuals versus the quantity 


log,,A.. Although the linear trends of EUS-WUS stations are perceptible, a 
10°c 


great amount of scatter remains. The slopes of the line fits are prescribed 


by the quantity “b" and the distance between the lines is the regional 


magnitude bias as prescribed by the value of "c'", The only flaw in the 


analysis, the assumption of fixed regional values of t*, does not seriously 


weaken our argument. We simply show the correlation between the regional 


magnitude bias and regional t* averages after allowing for approximate 


corrections for crustal effects. "C's" value and its confidence limits in- 


dicate that the numerical value of the reciprocal Am, (in magnitude units) 


cannot be less than At* (in seconds). 


Haskell, N. A., 1962, Crustal reflection of plane P and SV waves; Journal 
Geophys. Res., 67(12), 4751-4767. 
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Multipathing and local focusing have a significant impact on seismic 
wave amplitudes. Also, amplitudes of seismic waves at LASA (Mack, 1969) and 
NORSAR (Blandford, 1974; Chang and von Seggern, 1977) exhibit large systematic 
distance and azimuth dependent variations that can amount to a full order of 
magnitude, a figure that is probably true for almost any location on earth. 
At LASA these anomalies tend to cancel out if the data are evenly weighted in 
the azimuth-distance sense (Chang and von Seggern, 1977), but this cannot be 
assumed in general. Such amplitude variations may stem from multipathing and 
ray focusing and they cannot be predicted at individual sites without detailed 


structural knowledge and concomitant analysis such as 3-dimensional ray tracing. 


On the other hand, direct t* from slopes of spectral ratios appears to be 
considerably more consistent. Calculations of the crustal effect over a large 
variety of crustal structures on t* show scalloping of the spectra. If slopes 
are taken over a sufficiently wide frequency range, the t* measurement is not 
affected by most crustal structures because the overall slope of the spectrum 
does not change. However, individual structures could still change the slope 


somewhat. 


Studies of slopes in spectral ratios over many sensors at NORSAR (Der et 
al., 1977b) showed that while amplitudes of P waves varied by a factor of five 
giving a om, of several tenths, the corresponding standard deviation of At* 
was only .06 across the array. Thus, on the basis of relative accuracy alone, 
At* is a considerably more stable quantity than dm, - Plots of the standard 
deviations for single dm, and At* measurements as functions of mutual distance 
between station pairs further emphasize this point (Figures 3 and 22). While 
the standard deviation of dm, increases by a factor of 4.5 over the distance 
range the stations covered, the corresponding increase in the standard deviation 


of At* is by only a factor of 2.25. Also, t* starts at lower numerical values. 
Mack, H., 1969, Nature of short-period P-wave signal variations at LASA; J. 
Geophys. Res., 74, 3161-3170. 


Blandford, R., 1974, Short period signal-to-noise ratio at NORSAR; SDAC-TR- 
74-13, Teledyne Geotech, Alexandria, Virginia. 


Chang, A. C. and D. H. von Seggern, 1977, Study of amplitude anomaly and 


bias at LASA subarrays; SDAC-TR-77-11, Teledyne Geotech, Alexandria, 
Virginia. 
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That "b" is significantly different from zero and close to unity indicates 

that the crustal corrections are realistic and that log A, is a good overall 
estimator of the crustal effect. The remaining scatter (0 ~ .1 magnitude 
units) stems from: 1) uncertainties and possible biases of Bn, due to multipath- 
ing, focusing the errors in admittedly rough crustal corrections, and 2) the 
likely deviation of the individual station t* from the assumed regional values. 
Still, the terms associated with log) oA. and t* together account for about 74% 
of the total variance of Am, which, considering the inconsistencies of short 
period data, is a reasonably good result. The t* differential accounts for 53% 
of the total variance; the crustal effect for only 10%. Note that the EUS and 
WUS populations are almost completely separated in the Am, - log A, plane. 

The only overlapping point is MOID, a station located in a silo sunk into un- 
consolidated silt. 

The most critical parameters in crustal response calculations are the 
elastic constants and the density at the surface. Figure 24 shows a plct of 
1og1 9A, versus the quantity (1/2) Log, )(o,4,°), called the surface impedence 
term (6, and a, are the density and compressional velocity at the free 
surface). The correlation is quite good, suggesting that (1/2) Log, 9 (90,7) 
could be used instead of log,gA,- Figure 25 shows a plot of (1/2) 10g, 9(0 4, °) 
versus am» which again shows the separation of the EUS and WUS populations, 
while the distance between the two Straight line fits indicates about the same 


regional magnitude bias. In this case the t* coefficient was -1.33. 


The data set of Booth et al. (1974), without individual t* for each 
station, cannot measure the noise term ¢ in equation (2). However, it 
establishes jointly with our log, A, term that crustal effects are signifi- 
cant and that the EUS-WUS bias is not an artifact caused by crustal effects. 
The values of the coefficients "b" and "c" are consistent, on theoretical 


grounds, with the expectations. 


The pulse from a 50 kt explosion has teleseismic P waves similar in fre- 
quency content to teleseismic P waves from earthquakes: hence, b ~ 1. The 
Value of c = -1.33 can also be justified by synthetic seismogram simulations. 
Synthetic seismograms were computed for nuclear explosions of various yields 


in granite and then passed through causal constant~Q filters with various t*. 
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The results are shown in Figure 26. A line with slope -1.35 is also drawn for 
comparison. The slopes derived from the synthetics seem to be roughly equal 


to or greater than -1.35, depending upon the absolute t*. Another explanation 


for this result is that the empirical formula is derived from earthquake measure- 


ments while explosions have more high frequencies in their spectra. 


Another argument can be made for a single frequency component. The 


exponent in the attenuation law exp (-ft*) leads to a formula 


Am, ~ o-nfAt* 


1n10 (3) 
which yields 
Am, = -1.36 At* for f = 1 
and 
Am, =- 1.96 At* for f = 1.42 


(f = 1.42 is the average frequency of P waves at RKON). These values can be 
substantiated with spectral measurements. However, in the actual time 

domain magnitude calculation, correction for instrument response reduces 

the Am, This occurs because the dominant period is also changed by 
attenuation, and the instrument response changes rapidly as a function of 
period between 0.7 and 0.9 sec, which are typical P-wave periods at RKON 

and NTS, respectively. This effect is discussed in detail in the preliminary 


report of Phase I of the NTS experiment (Der et al., 1977). 


The NTIS experiment produced a set of At*, Am, and 1og, oA, values that 
are useful in testing the consistency of equation (2), the rms magnitude 
of the terme. The rms value of « based on the stations of the NTS experi- 
ment is quite large (o > .1 m.u.), an observation based upon the inconsistency 
(non-parallelisin) of corrected Am, and At* plots in Figure 4. The numerical 
values of the quantities plotted in Figure 4 are also summarized in Table II. 
Figure 27 is identical -o Figure 25, except that data from the NTS experi- 
ment have been superimposed on it. We used RKON as a reference point to 


bring the two data sets together. The figure shows that while FANV, GBNM, 
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and possibly SZNV, fit Booth et al's (1974) WUS population, all NTS points 

have higher my, values than the average WUS population. However, using t* 
values would be the way to determine anelastic energy loss for outgoing P waves 
from NTS. Still, because of the m results, we assume that there is focusing 
under NTS for events from the Northwest and Southeast where most of the events 
came from. By reciprocity (Chang and von Seggern, 1977), a network deployed at 
these sites would see higher magnitudes from NTS than from the rest of the WUS. 
We do not know what a differently deployed network would show. 
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NOISE STUDY 


To evaluate noise conditions at the SDCS sites we routinely measured 
the noise level preceding each event. The noise readings evenly covered 
each station's operational life, and the number of day and night time read- 
ings were roughly equal. Although somewhat biased by the seasonal distribu- 
tion of the operational life of each station, the gross distributions of 
noise amplitudes indicated the average noise conditions at each site. The 
cumulative distributions of noise log-amplitude readings for each station 
are given in Figures 28 through 40. The r2an values and standard deviations 


of noise are given in Table VI. 


Cumulative distributions were also computed for all available seasons 
for each station. These are shown in Appendix D. Table VI summarizes the 


means and standard deviations. Table VII lists station cocrdinates. 


Each station was also examined for daily noise variation. Figures 
41 through 47 show some selected plots of daily noise variations over 
several days. Much of this noise is obviously man-made, reflecting daily 
routine activity; generally, this noise is high frequency. One inter- 
esting example of this daily variation was at TQMS (Figure 41). the site 
of the SALMON nuclear explosion. There, the surface instrument showed a high 
daily variation of high frequency noise at the surface, but greatly reduced 
noise at the deephole instrument. The high frequency surface noise tended 
to follow the daily cycle of human activity. In quiet periods when “cultural 
noise" was low, the noise consisted mostly of common storm microseisms that 


had the same amplitude in tne borehole as on the surface. 


Daily variation at YFNV (Figure 43) shows great m 3-day maxima that 
must be associated with drilling activity at Yucca Flats. Some secondary 


peaks at midnight must also be associated with human activity, but we do 


not know the source of the noise. The sites FANV — OB2NV (Figures 42 and 44) 
also show peaks during the middle of the day that are presumably associnted | 
with normal working hours. Figures 45 through 47, that show stations 


GBNM, RXON and HNME, do not display the regularity in the daily noise 
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Cumulative freauency histogram of a11 noise readings at YF3NV. 


Figure 36. 
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a Geographical Coordinates of Stations Discussed in this Report 
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TABLE VII 


LAT 
46° 09" 43" 
40° 40° 20" 
272.13" 31" 
S713" By" 
37° 04' 06" 
37° 04" io" 
37° 04° 22" 
37° 04" 29" 
38° 38° 26" 


36° 41' 13" 
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LON 

57° 59! 

93° 40' 

116° 03! 
116° 03' 
116° 00' 
116° 00' 
116° 01' 
116° 02' 
116° 13' 


107° 13' 


09" 
20" 
28" 
28" 
07"' 
44" 
27" 
12" 
2a 


34" 
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variations. For GBNM, two peaks were picked during the days around 
Christmas; however, the peaks were not consistent with what was considered 


normal activity. 


The mean values of noise (in millimicrons) and logarithms and standard 
deviations of the logarithms are self~explanatcry. Seasonal variation is 
apparent at HNME, a station dominated by long-period storm microseisms. 
Stations at NTS show a similar seasonal variation. TQMS and GBNM on the 
other hand, show lower amplitudes in the Fall, which can be explained by 


seasonal variation of cultural noise. 
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CONCLUSIONS 


The results of spectral analyses show that when compared to RKON and 
HNME, a consistent loss of high frequency content in P waves can be 


observed at the SDCS staticns located in the WUS. 


‘ne average teleseismic P wave amplitudes are alse lower in the WUS, 

if crustal amplification corrections have been inade. Nevertheless, 
amplitudes of teleseismic P waves at NTS seem to be higher than at sites 
in other areas of the WUS. This is presumably due to focusing under 

the site. 

No indication of inherent bias in the determination of relative 
magnitude differentials was found. 


Relative t* derived from spectral ratio measurements is a numerically 
more stable quantity than the magnitude residuals. It has less scatter 
and depends less on interstation distance and the local crust. 
mend that local mantle magnitude attenuation is estimated from the 
empirical formula 


Am, = -1.36 At* 
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SUGGESTIONS FOR FURTHER WORK 


While one shortcoming of this study is the relatively small number of 
stations where both Mh and t* biases are known, the existing data set could 
be easily extended by adding more stations using historical LRSM daca. 

The multiple regression analysis involving mh t* and crustal amplification 
could then be used to better evaluate interrelationships between these 
quantities. Another useful project might include extending the study to 


SRO's and overseas observatories. 
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APPENDIX C 


Histograms and table of At* values for selected station pairs. 
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APPENDIX D 


Cumulative frequency histograms of noise readings at SDSC 
stations during each season. 
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